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The standard ( "fine-grained" ) interpretation of quantum density functional theory, in which den- 
sities are specified with infinitely-fine spatial resolution, is mathematically unruly. Here, a coarse- 
grained version of DFT, featuring limited spatial resolution, and its relation to the fine-grained 
theory in the L 1 n L A formulation of Lieb, is studied, with the object of showing it to be not 
only mathematically well-behaved, but consonant with the spirit of DFT, practically (computa- 
tionally) adequate and sufficiently close to the standard interpretation as to accurately reflect its 
non-pathological properties. The coarse-grained interpretation is shown to be a good model of formal 
DFT in the sense that: all densities are (ensemble)- V-representable; the intrinsic energy functional 
F is a continuous function of the density and the representing external potential is the (directional) 
functional derivative of the intrinsic energy. Also, the representing potential v[p] is quasi-continuous, 
in that v[p]p is continuous as a function of p. The limit of coarse-graining scale going to zero is 
i-C . studied to see if convergence to the non-pathological aspects of the fine-grained theory is adequate 

to justify regarding coarse-graining as a good approximation. Suitable limiting behaviors or intrinsic 
energy, densities and representing potentials are found. Intrinsic energy converges monotonically, 
coarse-grained densities converge uniformly strongly to their low-intrinsic-energy fine-grainings, and 
L ' + L°° representability of a density is equivalent to the existence of a convergent sequence of 
coarse-grained potential/ground-state density pairs. 
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Underlying electronic density functional theory (DFT) [THE] in both the dominant Kohn- 
Sliani <> and Orbital-FreefJi] variants, is a density functional, F[p], representing the mini- 
mum kinetic plus (Coulomb) interaction energy of N electrons compatible with the density 
p. Clearly, the properties of F are of great importance, and the formal development of DFT 
seems, at least implicitly, to involve assumptions that (a) the intrinsic energy functional 
is differentiable, (6) each density can be selected as a ground state density of an external 
potential, which is the functional derivative of F at p, (c) F is continuous. None of these is 
true. Better to say: none of them holds in the standard interpretation. For, there are things 
to be interpreted. Continuity requires a topology, and derivatives come in various types. 
Might the interpretational task extend even to the terms 'density' and 'potential'? The 
standard interpretation is fine-grained in that a density is assigned to every point ( "almost- 
every" point, technically). If we construct a regular partition of space into cells and regard 
number-in-cell divided by cell- volume as 'density', then we arrive at a coarse-grained, or 
resolution-limited, interpretation. The thesis of this paper is that this coarse-grained model 
is consonant with the spirit of DFT, mathematically benign, and a good approximation 
to the fine-grained theory in that it faithfully reflects its non-pathological aspects. The 
coarse-grained interpretation renders every density the ground state of some potential, es- 
sentially uniquely (Hohenberg-Kohn theorem), it satisfies assumptions (a-c) above, and the 
representing potential v[p] is 'quasi-continuous' in that v[p]p is L 1 continuous. The good 
approximation properties are phrased in terms of limits as the coarse-graining scale is taken 
to zero. Intrinsic energy converges monotonically to its fine-grained value. A fine-grained 
density is V-representable by a potential in L 3 / 2 + L°° if and only if there is a sequence of 
coarse-grained densities which converge in an appropriate way, along with their representing 
potentials, to the fine-grained data. Coarse-grained densities approximate in a uniform way 
all the fine-grained densities of low intrinsic energy which project to them. 

To expose a little better the basic idea and its consonance with DFT, consider the 
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constrained-search formulation [1JJ-[14|]. The focus of our interest is the Valone-Lieb intrinsic 
energy functional J12|, ll3( defined as 

F[p] :=inf{Trr(T + K e ) : T ^ p} . (1) 

Minimization is over antisymmetric mixed states of N identical fermions T = ^\ 
with Xi = 1 having single-particle density p (this relation is denoted T p). T is the 
kinetic energy and V ee the interaction energy among the fermions, so that TrT(T + V ee ) is 
the expectation value of kinetic-plus-interaction energy in the state T. We have in mind 
electrons in particular, and the interaction is assumed no more singular than Coulomb. 
States of the system are partitioned into equivalence classes according to their density. 
Ultimately, only the lowest intrinsic energy state in each equivalence class matters (assuming 
that the minimizer exists). For external single-particle potentials v in some set V, we define 
the ground-state energy as a Legendre-Fenchel transform of F: 

E[v] := inf |f[p] + J v(x)p(x)dx : p e X>| . (2) 

If T + V ee +v really has a ground state with a density in T>, this will pick out its energy. The 
point is that the potential couples to the density only, so that the minimization problem 
"find the ground state" can be split: first find the minimizing p in Eq. @, then go back 
and identify the state from Eq. ([I}. 

Suppose now a partition of space into disjoint cells, and take V to consist of potentials 
uniform over each cell CP-measurable) . In that case, the potentials are sensitive to the 
average density over a cell, but not to the variation within; that variation still exists, but 
we relinquish control over it. Densities fall into equivalence classes, coarse-grained densities, 
identified by average density over each cell. Ultimately, only densities (possibly nonunique) 
of lowest intrinsic energy in each equivalence class matter. With V thus adapted to *p, we 
make a reinterpretation, taking T> to consist of coarse-grained densities, and Eqs. ([T][2]) apply 
as they stand. This is the basic idea of coarse-graining. We are not restricting the densities 
or the states or changing the quantum mechanics by working on a lattice, for example, but 
are working with a limited spatial resolution, and thus cannot distinguish all densities. It 
is in the nature of the problem that the fine-grained member of each coarse-grained density 
which is selected is the one with lowest intrinsic energy. The projective character of the 
coarse-graining on the density side of the ledger is naturally paired with a restriction on 
the potentials. The partition *p has logical priority. It determines a collection of number 
observables which in turn determine the denotations of 'density' and 'potential'. 

Some densities are ground-state densities in some external potential (V-representable) . If 
v is one such potential, adding an arbitrary constant to v preserves this relationship. The 
Hohenberg-Kohn theorem (which holds also in the coarse-grained theory, see Section [3]) says 
that this is the only freedom. We establish a convention for fixing the constant: v[p] denotes 
the potential having p as a ground-state density with the constant adjusted so that 

F[p] + (v,p) := F[p] + J v(x)p(x)dx = 0. (3) 

Equivalently, = 0. This particular convention is somewhat arbitrary, but some such 

in necessary for all the statements in the introductory paragraph to make sense. 

The body of the paper unfolds as follows. Section [2] discusses relevant aspects of the 
fine-grained, theory. In particular, the "bad-behavior" alluded to earlier, which is not so 
well-known, is exposed. This provides both motivation and context for the coarse-grained 
model. A conclusion is that the difficulties are connected to short distance scales, which 
suggests that some sort of short-distance regularization is called for. As Chayes, Chayes 
and Ruskai 15] put it, "long-distance difficulties do no t g enerally occur in HK theory." But 
short-distance ones certainly do. Lattice models |15l - fl7j are one means of short-distance 
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regularization. But such discretizations, if taken seriously, involve fundamental changes to 
quantum mechanics. Their relation to the continuum theory, and the continuum limit, is 
tricky and ambiguous. The projective character of the regularization achieved by coarse- 
graining is much gentler. One might say that the modification is at an epistemic rather 
than ontic level. It is for that reason that we use the somewhat awkward term 'fine-grained' 
rather than 'continuum'. The continuum nature of space is recognized by the coarse-grained 
interpretation. Some technical results are also developed for use in Section [5J 



As for lattice formulations, approximations which use a limited basis for states or density 
matrices (studied in great detail in the series of papers [l8l - [2lj ]) must be differentiated from 
the current approach. A coarse-grained density represents an equivalence class of infinitely 
many densities. Insofar as low energies are of interest, most of them are not very relevant, 
but the selection is an energetic one, not an a priori choice. Coupled with the projective 
treatment of densities in the coarse-grained model is an injective approximation of densities; 
it is here that a functional palette is limited. This might seem a somewhat uncomfortable 
aspect. We would like to handle Coulomb potentials of point charges, for instance. But 
an atomic nucleus is not a point charge and if we treat it as such, we are half a step from 
coarse-graining at the femtometer scale, anyway. 

The coarse-grained model is set up in Section [3] and results from a previous paper[22] 
are reviewed. The new results for the single-scale coarse-grained model are discussed and 
proven in Section 21 showing that it is free of the bad behavior which afflicts the fine-grained 
theory. Theorem 14.11 shows essential continuity of the intrinsic energy, Theorem 14.21 shows 
that physical directional derivatives of F[p] coincide with representing potentials v[p] and 
Theorem 14.31 shows that v[p]p is continuous. Kohn-Sham theory is taken up in Section 
[5j where the exchange-correlation potential is shown to exist at least as the directional 
derivative of the exchange-correlation energy. 

Sections [3] through [5J show that the single-scale coarse-grained model is a well-behaved 
model of formal DFT. Its acceptability, however, depends to some extent on its ability to 
reasonably approximate the fine-grained theory as the coarse-graining scale goes to zero. 
This is the subject of Section [51 where a multi-scale coarse-grained model is set up. Taking 
the coarse-graining scale to zero is naturally thought of as a process carried out in time, 
and the multi-scale model is ultimately little other than this process viewed sub specie 
aeternitatis. It is a convenient tool however, and serves to remove the impression that a 
fundamental length scale is inherent in the coarse-grained interpretation. A basic result is 
that all the intrinsic energy of a fine-grained density is recovered monotonically in the limit. 
Representability of a density by a L 3 ^ 2 + L°° potential is reliably and faithfully signalled by 
the coarse-grained model. Proximity of coarse-grained densities to the low-intrinsic-energy 
fine-grained densities in the equivalence classes they represent is also examined. Among 
other things, these results are taken as validating the claim that the coarse-grained model 
approximates what it is supposed to. 



A nodding familiarity with basic Banach space functional analysis is presumed in the body 
of the paper. Appendix lAlcontams a brief review of concepts, definitions and notations which 
may be helpful to readers needing a quick reminder, or more. 



A couple of delimiting remarks are in order before we begin. In this paper, mixed states are 
always allowed. So, all statements about V-representability refer specifically to ensemble-V- 
representability. The results are applicable to abelian spin-density-functional theory, where 
up-spin and down-spin are roughly separate species. However, genuine non-abelian spin- 
density-functional theory. fTA [23 . [24 [ treating all components of the spin, is unfortunately 
well beyond the scope of this work. 
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2. FINE-GRAINED DFT AND ITS DISCONTENTS 



This section is an essay on fine-grained DFT, focussing on the difficulties mentioned in 
the Introduction. Lieb's landmark 1983 paper[13[ formulates the theory in the precise form 
in which we will consider it. The recent article 25] of van Leeuwen and the bookQ by 
Eschrig are also recommended, though some assertions found in them are here rejected. 
Further background on convex and nonsmooth analysis can be sought in the books [261429] ] . 
Chapter 5 of Aubin and Ekeland[27] contains an interesting and idiosyncratic discussion of 
the Ekeland variational principle. Appendix [A] contains a refresher on functional analysis 
and Lebesgue spaces which is relevant. 

It is demonstrated in §1 of Lieb's paper[l3[ that the effective domain of F is 

J N := \p : p(x) > 0, p 1/2 E H\R) 3 , [ p(x) dx = n) . (4) 



Here, fl" 1 (R 3 ) is the Sobolev space of functions / such that J |/| 2 + |V/(a;)| 2 dx < oo. The 
gradient can be interpreted in terms of Fourier transform; classical differentiability is not 
a requirement. The remarkable fact which makes Jjss the right domain is that there is an 
A-particle wavefunction tp with finite kinetic energy and single-particle density p if and only 
if p G .7 m [l3l [30l - [32l | . "Finite kinetic energy" is understood in quadratic form sense, i.e., 

Si!=i I l^iV'l 2 < °°j so V 2, does not need to be square-integrable, just V-0. This condition 
is sufficient to guarantee that also the Coulomb interaction energy is finite. Some pertur- 
bation, for example by an arbitrary bounded interaction, would not affect conclusions. We 
assume without loss that the interaction is bounded below by zero. Since J7at is convex [l3j. 
forming mixtures will not produce densities outside Jn, so it is the right domain whether 
considering mixed states or pure states. (A set X is convex if, whenever x and y are in X, 
so is the line segment Xx + (1 — X)y, < A < 1.) 

F[p] is now defined as in Eq. ([T]), with V — Jn- F has two important properties which 
allow the analysis to move forward: 

a. F is convex: F[sp + (1 - s)p'] < sF[p] + (1 - s)F[p'} for < s < 1. 

b. F is lower semicontinuous with respect to the L 1 topology. 

Assertions about ways in which nearby densities are similar requires a well-defined notion 
of 'nearby,' i.e., a topology (see Appendix [A"| . Different topologies can be appropriate and 
useful in different ways. The statement about lower semicontinuity means: Given p € 
and e > 0, there is d > such that \\p' — p\\i < d implies that F[p'] > F[p] — e. The name 
should now be clear. // the conclusion had instead been F[p'] < F[p] + e, that would be 
upper semicontinuity of F. Continuity is the conjuction of lower and upper semicontinuity. 
Why is F not upper semicontinuous? That is connected with the important inequalies[13| 

c\p 1/2 \ 2 m<F[p]<c' + c'V /2 \ 2 m, (5) 
where the squared Sobolev seminorm is 

a/212 _ _ 



\P 



\ 2 m = J\Vp^\ 2 dx = J^dx, 



which will also be recognized as the von Weiszacker term[2j. The constants in the upper 
bound can depend on the interaction, but the lower bound depends only on the kinetic 
energy. It is a precise form of the physical intuition that strong oscillations in the density cost 
kinetic energy. Consider, for instance the density p'(x) ~ [1 + er/(x) sinx/£] 2 p(x), where 

< 77 < 1 is a smooth function with bounded support. \\p' — p\\i — 0(e), but \p' |^i — 
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0(e 2 /£ 2 ) as I — > 0. This is the first serious short-distance difficulty. Such problems are 
a major motivation for the development of the coarse-grained approach. Although this 
discussion was framed with respect to L 1 for concreteness, it is easy to see that it also 
applies to L 1 n W for 1 < p < 3. 

Now we embed Jn in a Banach space A with a norm at least as strong as L 1 (meaning 
the norm dominates some fixed multiple of the L -norm). Lieb noticed that a Sobolev 
inequality 

IHIs^cd^l^ + yii) (6) 

combines with (J5J) to show that Jn is actually contained in L 3 as well as L 1 . Consequently, 
he chose X = L 1 n L 3 with norm || • \\x = \\ • ||i + || • 1 1 3 . But there is freedom here. 
Any L 1 n L p with 1 < p < 3 would work, as would L 1 . However, a norm at least as 
strong as L 1 is important to maintain lower semicontinuity of F. The move to a Banach 
space puts us in position to use some relevant ideas and results of convex analaysis, and 
is suggested by simple observation that an external potential v acts as a linear functional 
(v, p) :— J v(x)p(x) dx. F is extended by taking F = +00 off J^. This preserves convexity 
and lower semicontinuity [3]. The subset of X on which F < +00 is called the effective 
domain of F and denoted dom F. Of course, dom F is simply Jn , but we will use the former 
notation sometimes because the effective domain in the coarse-grained model is different, 
but plays the same role. Note that J7}v is contained in the closed affine hyperplane 

X N : = |p G X : J p(x) dx = AT J . 

Since elements of X, and even of Ajv, can be negative, we refer to them generally as 'quasi- 
densities,' reserving 'density' for the non-negative elements. 

The ground-state energy E[v) is defined according to Eq. @. But, for which potentials? 
The expression on the right-hand side of Eq. ([2]) makes sense for many potentials, and we 
will return to that point. But insofar as we want F[p] and E[v] to form a Legendre-Fenchcl 
transform pair, we should concentrate on its restriction to the dual space A* of A. E, as a 
function on A* is concave and upper semicontinuous from its definition, and general results 
of convex analysis (due to the known convexity and lower semicontinuity of F), guarantee 
that 

F[p]= sup {E[v]-(v,p) : v€X*}. (7) 

There is an apparent asymmetry here due to our choices of signs, but also a more glaring 
asymmetry: L 1 (M 3 ) is not the topological dual space of L°°(M). This is repaired by using 
the weak and weak-* topologies on A and A* (see Appendix [5} under which they are 
topological duals of each other. Symmetry is thus restored. Just as important, F and E 
are lower and upper semicontinuous respectively with respect to these topologies due to the 
remarkable fact that for a convex subset of a Banach space, closure with respect to the norm 
or weak topology are the same. 

Lieb worked with A = L 1 D L 3 . That has the advantage over L 1 that the Coulomb 
potential of a point charge is in the dual space A* = L 3 / 2 + L°°, but it seems we need 
to go outside A* to find all the potentials which may be of interest anyway. Consider, for 
example, the harmonic oscillator potential \x\ 2 . It is not in any of our A*s, yet it does 
not really pose any particular difficulty. Eq. @ works for it, as long as we restrict the 
minimization to p G Jjq . [Outside that set, the right-hand side of Eq. @ involves nonsense 
such as co + (—00).] It begins to appear already that y-representability by a potential in 
A* has a special status. We will call this restricted notion X* -representability. 

If we suppose that p is V-representable, that is, a ground state density for some potential 
v[p], and that F[-] + {v, •) is somehow smooth, Eq. ^ suggests that the representing 
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FIG. 1. Varieties of local behavior of F[p] along linear slices in Xn, through a density p. In (b) 
and (c), F takes the value +00. 



potential is some kind of functional derivative of F: SF/Sp = — v[p\. Smoothness is not in 
the cards, but convex analysis offers the derivative-like notion of subdifferential that does 
not depend on it. The subdifferential of F at p, traditionally denoted dF, is a subset of X* , 
the elements of which are called subgradients. The definition is 

v G dF[p] <=> F[p'\ > F[p] + (v,p' - p), Vp' G X. 

A subgradient is the 'slope' of a continous hyperplane which touches the graph of F at p, 
but is nowhere above it. This is a looser notion than the ordinary one of tangency: the 
subdifferential of an ice-cream cone at its point has many elements. The relation @ shows 
that dF(p) consists precisely of those potentials for which p is a ground-state density. From 
the Hohenberg-Kohn theorem, we deduce that the subdifferential at p either is empty, or 
its elements differ only by a constant. 

How does the subdifferential relate to other notions of derivative? Probably the simplest 
such is that of directional derivative. The directional derivative of F at p G domF in the 
direction dp £ X is 

F'{p; 5p] = lim - (f[ P + s Sp] - F[p}) , (8) 

if the limit exists. But, convexity guarantees that the difference quotients are nonincreas- 
ing as s — ¥ 0, so the limit always exists, though it may be infinite. In particular, if 
F[p + sSp] = +00 for all s > 0, then F'[p;Sp] — +00. The domain of the directional 
derivative at p, denoted domF'[p; •], is defined as all 5p such that F'[p;Sp] < +00 If, for 
every Sp, the directional derivative is finite and satisfies F'[p;—Sp] = —F'[p;Sp], then the 
various directional derivatives fit together into a linear functional. If this linear functional 
is continuous, it is called a Gateaux derivative, or G-derivative. (Beware, the continuity 
requirement is not always imposed 29j.) If p is X*-representable, the variational principle 
assures us that 

F'[p;Sp] >- J vSpdx. 

An important question then is to what extent equality holds. There are certainly 5p for 
which F' [p; ■ ] is infinite. 

The situation is illuminated by looking at varieties of local behavior of F along line 
segments in X. They are depicted in Fig. [T] First, of course if / 5pdx ^ 0, then for any 
s 7^ 0, p+sSp is outside Xn and F[p+sSp] = +00, so this is a trivial case of Fig. Q](c). But, 
even if the perturbation does not stray from Xn, (b) and (c) occur for every p G domF. 
This may be surprising at first. But it is not so hard to find two sources of this kind of 
behavior. First, there are perturbations for which p + s Sp, for s on one or both sides of 
zero, falls outside X^, the set of elements of Xn which are everywhere non- negative. For 
example, this will happen for s > if the negative part of Sp falls off too slowly at infinity, 
depending on p. If the positive part also falls off too slowly, then we have a picture like (c). 
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Secondly, there are perturbations for which p+s Sp is in for s in some open neighborhood 
of zero, but not in Jn- This, too, leads to case (c). A severe enough discontinuity of Sp, 
for example, will do the trick. A variation on the construction we used to show that F is 
locally unbounded on Jn provides a more subtle example. Adding oscillations of the right 
amplitude to p at a hierarchy of ever shorter wavelengths will produce a sequence converging 
with respect to A-norm to p + 8p G Xi, but with an infinite kinetic energy. 

Falling outside X^ under perturbation produces both scenarios (b) and (c). Falling 
outside Jn while remaining in X^ gives only scenario (c). To see this, consider p, p' £ Jn, 
and suppose ||p'/p||oo = M > 1 is finite. Then, we already know that p\ = (1 — X)p + \p' 
is in for < A < 1. And, p\ is in X~^ for —M < A < 0. For A in this latter range, p\ is 
in Jn if \^/p\\h 1 is finite. Since 

^ = ^((1 - A) 2 |VH 2 + A 2 |Vpf + 2A(1 - A)Vp • Vp' 
Ap x Ap\ V 

and |2A(1 - A)Vp • Vp'\ < |A(1 - A)|(|Vp| 2 + |Vp'| 2 ), taking A = -1/(2M) gives p : = 
P-1/2M > p/2 > p'/{2M), so that 

|Vp| 2 lVp| 2 nnT \V P '\ 2 
—rz- < Q-^r- + 2M 1 „ , < oo. 

But, — 1/(2M) is halfway to the critical value of A, and nothing prevents taking that step 
again and again. Thus, if p and p' are in J7jv, the interior of the line segment {(1 — X)p+ Xp' : 
AsKJfl X^ is also in Jjv ■ This argument does not tell us what happens at an end-point 
of this segment. However, lower semicontinuity of F implies that either it is in J7}v , or F 
diverges to +00 there. 

This second class of density perturbations is apparently another short-distance problem, 
and will be absent in the coarse-grained model. But, the first kind, connected with the fact 
that every density has a tail which dies away, will still occur there. On a slice entirely in 
domi* 1 , we will get something like Fig. Q] (a). We would hope that there was no kink, but 
there is no proof that is the case, even if F has a (necessarily essentially unique) subgradient 
at p. 

Here is a tempting argument put forward some time ago[33|, which unfortunately keeps 
recurring. It purports to show that at an A*-representable density, a conventional Gateaux 
functional derivative exists and coincides with the representing potential. The epigraph of 
F, the set of points on or above its graph, 

epiF = {{p,z) e X x K : F[p] < z} , 

is convex. Any directional derivative gives a line disjoint from the interior of epi F, hence 
can be extended to a full hyperplane by the Hahn-Banach theorem (the interior is the set 
of points which have a full open neighborhood contained in epi-F). Since such a hyperplane 
must coincide with the subgradient by assumption, the conclusion follows. This does not 
work. We have just seen that with Fig. [T] A Gateaux derivative does not give infinite 
directional derivatives. That same figure also shows the flaw in the argument: the interior 
of epi F relative to Xn is empty. 

Of course, the only directional derivatives F[p; Sp] which have a real physical significance 
are those for which both p and p + s Sp for some s > are in dom F. The others are an 
artifact of the Banach space imbedding, and are sure to be ill-behaved, so that the best that 
could be expected is 'quasi-G-diffcrcntiability,' meaning that all those physically significant 
directional derivatives fit together into a continuous linear functional. 

Now we turn to the 1^-representability issue. Which densities are ground state densities 
of some potential (y-representable)? A coherent answer to the question probably requires 
deciding what will count as a potential. One possibility, motivated by the Legendre-Fenchel 
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duality, is to consider only v's belonging to X*. This X*-representability, already intro- 
duced, is a somewhat restricted notion but it is well-structured and something can be said 
about it based on corollories of the celebrated Ekeland variational principle [jyj. These 
results will be used in Section [6] From Eq. ([2]) we have 

E[v]<F\p] + (v,p), v€X*,peX. 

If p and v saturate the inequality, then v £ dF(p). Suppose that it is almost saturated, so 

F[p] + (v,p) <E[v]+e. (9) 

In that case, we say that v is in the e-subdifferential d e F[p\. Physically, this means that p 
almost achieves the ground-state energy in v. Does that imply that it is almost a ground- 
state density of v? Not quite, but the Ekeland variational principle does warrant (for 
example, Thm. 1.6.2 in Ref. [26|) the following: For any A > 0, there is a density-potential 
pair p\ € domF and v\ £ X* , such that p\ is a ground-state density of v\ and 

F[px] + {v,px) <F[p] + {v,p), 
\\p\ - p\\x < A, 

\\vx-v\\ x * <e/X. (10) 

(Note that v is on both sides of the first inequality.) Interesting conclusions can be drawn 
from this. Given a density p, pick any potential and find e to satisfy inequality ©■ Then, 
taking A = 1/n, the theorem gives us a sequence of X*-representable densities converging 
to p, showing that X*-representable densities are dense in domF. This is the Br0nsted- 
Rockafellar theorem. But, note that control over the norms of the representing potentials 
gets progressively worse as n — > oo. Another interesting choice is to set A = -Je. In that 
case, \\px ™ p\\x < an d \\vx — v\\x* < so that if e is small only small perturbations of 
both p and v to find a ground-state density /representing potential pair. This does not say 
that all nearby JT*-representable densities have nearly the same representing potential. 

In thinking about the relationship between X*-representable densities and their associated 
potentials, it is useful to work in the product space X x X* . The graph of the subdifferential, 

graph := {{p, v) e X x X* : v £ dF(p)} , 

then naturally suggests itself as an object of study. If a pair (p,v) is not in graph dF, then 
p is not a ground state density for v. This is a little crude. It might be even better to have 
a function which measured by how much it fails. To this end, we introduce the energetic 
excess of p with respect to v, 

A[p,v]:=F[p] + {v,p)-E[v}. (11) 

A[p, v] is the lowest energy attainable with quasi-density p in the presence of v, relative to the 
ground state energy of v, and therefore A : X x X* — > [0, oo]. Clearly, graphed = A _1 (0). 
Note that, although the natural description of the significance of A seems asymmetric, p 
and v are really on essentially equivalent footing in A[p, v] as follows from the conjugacy of 
F and E. 

Some properties of A follow immediately from those of F and E which have already been 
discussed. Namely, it is convex separately in each of p and v since the pairing (v, p) is linear, 
and it is lower semicontinuous with respect to either (weak) x (norm) or (norm) x (weak-*) 
topology on X x X* since (v, p) (v, p) is continuous with respect to either of these. 
This means in particular that graph <9F = A _1 (0) is closed with respect to either of these 
topologies. In fact, if only 



A[p n ,v n ] -> and (p n ,v n ) -t (p,v), 
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then 

(p,v) G A _1 (0), F[p] = lim F[p n ], and £?[«]= lim E[v n ]. 

n— >oo n— >oo 

The last two limits follow since -F can only decrease at the limit and E can only increase, 
but A > 0. 

The bonus convergence of F and E implies that the subset of graph dF which respects 
the gauge convention of Eq. (J3J is itself closed. And, it further follows that the map p H> 
||u|jo]||x*) hke F[-\, is lower semicontinuous, where for convenience we declare ||i>[p]||x* = 
+00 if p is not X*-representable. For, suppose that v[p] exists but ||f["]||x* is n °t lower 
semicontinuous at p. Then, there is a sequence pi,p2, ■ • ■ in the domain of v[-] such that 
ll u L°.7']||jf* — ll 1 '!/ ]!!^ * — e - Since X* is a Banach space dual, bounded subsets are weak-* 
compact, which means that there is a subsequence w[p nj ] which is weak-* convergent to vq 
with norm not exceeding ||u[jo]||x* — e - But that would mean that vq is the representing 
potential for p respecting the gauge-convention, which is clearly impossible since its norm 
is too small. 

Part of the conclusion drawn from Ekeland's variational principle [see (fTU|)] takes a curious 
form when written in terms of the energetic excess A. We equip X x X* with the norm 
|| (x, w)\\xxx* = \\ x \\x + II HI jf"*- Since this is stronger than the topologies discussed in the 
previous paragraph, the convergence results stated there still hold. Now, if A[p, v] < e, then 
there is (p',v') e A _1 (0) with \\(p',v') — (p,v)\\xxX* < Thus, the minimum of A is 

necessarily "narrow" in some sense. 

Beyond X*-representability, not much can be said in general because not much is known. 
Englisch and Englisch[35| produced some examples of densities (in Jm) which are not V- 
nonrepresentable by any potential which is a function. This shows at least that such densities 
exist. The examples are for a single particle and use the fact that the only candidate 
potential in that case is given by (V 2 y / p)/y / p. These examples were further analysed by 
Chayes, Chayes and Ruskaifl5j. One is p = (a + b\x\ a+1 / 2 ) 2 for x near zero with < b < a 
and < a < 1/2. Non-V-representable densities can be undeniably physical, as shown by 
the examples [4j of single-particle excited states in a central potential with a node in the 
radial wavefunction, p ~ (r — ro) 2 . The non-V-rcprcsentability of all of these examples is 
clearly due to short-distance problems, which does not prove lack of other sorts, but provides 
more motivation for short-distance regularization. 

Thus, X*-nonrepresentable and even V-nonrepresentable exist, but how common are 
they? Suppose we equip J7jv with the metric d(p,p') = \\p — p'\\x + \F[p] — F[p'}\. Prop. 
[C]m Appendix [Cl shows that p4 H^MHx* is nowhere upper semicontinuous, and this has 
the consequence that J \"*-nonrepresentable densities are topologically generic in J7jy with 
respect to the topology induced by the metric d. In this sense, most densities in J7jv fail 
to be JT*-representable. The demonstration of this involves potentials which oscillate on 
extremely short wavelengths, so this is again a short-distance issue. However, one should 
note that only low intrinsic-energy densities are involved (at least explicitly). 

It might be suggested that absence of V-representability over all of dom F is disquieting, 
but far from a disaster. Ultimately, we are only interested in the V-representable densities, 
whichever ones those might turn out to be. The problem becomes more acute in Kohn-Sham 
theory. The Kohn-Sham decomposition of the intrinsic energy is 

F[p}=T s [p]+E H [p]+E xc [p}. (12) 

Here, T s [p] denotes the intrinsic energy for a non-interacting system of quasi-density p. Note 
that this has the same domain as does F. The Hartree energy Eh[p] is simply the classical 
Coulomb energy of the charge distribution — ep. The exchange-correlation energy E xc [p] is 
then defined by the equation. The point of this decomposition is that it facilitates a sort 
of self-consistent field approach. One solves a problem for non-interacting particles with a 
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potential which is the sum of external, Hartree and exchange-correlation potentials, self- 
consistency being imposed on the latter two. This has proved a very successful strategy for 
practical computations. If p is a ground-state density of an interacting system in external 
potential v[p], then it should be a ground state density of the non- interacting system in 
external potential 

v s [p] = v[p] + <f>[p] + v xc [p], (13) 

with <fi[p\ — SEn/Sp and u^p] = 8E xc /5p (note the sign convention). This raises a cou- 
ple of problems. What if p is (interacting) V-representable but not non-interacting V- 
representable? Then, evidently, the solution does not exist. Consequently, it is very im- 
portant to the Kohn-Sham enterprise that the sets of V-representable, or maybe we should 
say X*-representable, densities for the interacting and non-interacting systems coincide. At 
a density in the intersection, Eq. (fT3|) can be used to define i^d/)]- Even then, however, 
the existence of a functional derivative of E xc and its coincidence with v xc is not assured. 
Subdifferentiability of F and T s is not strong enough for that. 

Since it is the bad behaviors of the fine-grained theory which interest us, let us summarize 
them. F is not continuous in V-norm, even restricted to its effective domain. There are 
densities that are not V-representable, and many more that are not A*-representable. This 
puts two blocks in the way of well-definedness of a functional derivative of F . First the non- 
universality of subdifferentiability, and second the lack of a demonstration that directional 
derivatives coincide with a sub-gradient. All of these seem at least partially connected to 
abnormal occurences at asymptotically short distance scales. This motivates the coarse- 
grained approach, wherein the short-distance scale degrees of freedom are allowed to relax 
energetically. 



3. COARSE-GRAINING: NOTATION AND BACKGROUND 

In this section, we turn our attention to the coarse-grained model. Much of it will formally 
resemble the fine-grained theory outlined in the previous section. The term 'density' is 
intended in the coarse-grained sense unless otherwise noted. 

As explained in the Introduction, a (coarse-grained) quasi-density is a collection of number 
observables p(R)vol(R), one for each cell R of a partition *p of R 3 , satisfying a summability 
condition (see below). For purposes of working in a vector space, the "number observables" 
are allowed to have any real value. Thus, a quasi-density is identified with an equivalence 
class of elements of X , where two elements of X are equivalent if they have equal integrals 
over every cell. Ultimately, the only members of such an equivalence class which are essential 
are those which belong to Jn, if there are any. Still, a useful and easily visualized surrogate 
for a quasi-density is the everywhere-defined function which is uniform over each cell R CP- 
measurable), and equal to p(R). This is a special element of the equivalence class associated 
with p as we have defined it, but it far too discontinuous to be in J^. Some variation in 
the sizes and shapes of the cells of is permissible, but it must be limited (see Ref. (HI)- 
The partition generated by any sort of regular lattice is acceptable, and for simplicity, we 
keep a cubical lattice in mind. 

For a real-valued function / on *p, define the L 1 norm ||/||i = X)jte<p l/(-^)l va K-R) = 
/ |/(t)| dr. The Banach space of such functions with finite norm is denoted X. This is the 
space of quasi- densities. Certain subsets of X are given names. The set of everywhere non- 
negative elements, that is, the (proper) densities, is denoted X + ; that of those everywhere 
strictly greater than zero is denoted X ++ . Subsets of elements which integrate to a particular 
value are indicated with a subscript, e.g., denotes the properly normalized densities, 
and legitimate density perturbations are located in Xq. It is the nature of the normalization 
condition in particular which suggests the L 1 norm. Although it seems fairly natural, one 



11 



might question its appropriateness. Most of the results actually depend only on the topology 
associated with this norm, and one should note that it is the weakest topology making the 
number-in-cell observables and the total number continuous. More discussion of this point 
can be found in Section f6 Al 

The intrinsic energy is defined as in Eq. ([1]). Since the equivalence class in X associated 
to p E X intersects Jn if and only if p £ X ++ , the effective domain of F is 

domF :={peX : F[p] < +00} = X+ . 

The infimum in the definition of F[p] is guaranteed to be attained [22^. but there is no 
guarantee that there is only one fine-grained density which attains it. In some rough sense, 
domF here is a much bigger subset of X than Jn is of X. As in the fine-grained theory, F 
is convex and lower semicontinuous [22j . One thing coarse-graining has achieved is a bound 
for F: 

< F[p] < iVF max , P eX+. (14) 

NF mlix is the energy required to pack all N particles into a single cell (maximum such if 
the cells are not identical). 

The main result needed from Ref. [22j concerns unique and universal V-representability. 
As discussed in the Introduction, potentials are constant on cells. X* , the dual space of X, 
consists of CP-mcasurablc functions which are uniformly bounded. But we must go outside 
that space to find many of the needed potentials. The interesting set, that of ^-measurable 
functions bounded below but not necessarily above, is denoted V. 

Theorem 3.1 (Coarse-grained Hohcnbcrg-Kohn). For p 6 X ++ , the infimum in the defini- 
tion of F[p] is attained (at mixed state 7, say) and there is precisely one potential v[p] G V 
with a ground state (j) of density p. 

Proof. See Ref. Q. □ 

The constant offset in the potential is fixed here by the convention of Eq. (J3)): 

(T + V ee +v) 7 =E[v[p}]=0. (15) 

This has the effect that all external potentials associated with densities are bounded below 
by -F max (see Eq. [M]). 

4. SINGLE-SCALE MODEL 

In this section, the new results on the single-scale coarse-grained theory are studied. In 
Section \4 Al the Theorems are stated, and discussed with some indication of the methods 
of proof. Detailed proofs are given in Section 14 Bl 

A. ideas 

It is cumbersome to have to maintain strict normalization of densities and states at all 
times, so we will work with a slight modification of F. F[p] is defined as in Eq. ([1]), but with 
the normalization restriction on states lifted. Thus, F agrees with F on X^ (which is the 
physically important set). But, for p G X + , it scales linearly with the normalization, -F[A/?] = 
A-F[p], and is +00 elsewhere. Note this maintains convexity and lower semicontinuity. The 
theorems are stated in terms of F, but their translations into terms of F are easy. 

As seen already, coarse-graining renders F bounded in a fairly trivial way. It is actually 
even continous with respect to the L 1 topology on its domain, X + . 
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Theorem 4.1. F is continuous on X + with respect to the L 1 topology. 

Since F[p] is already known to be lower semicontinuous, proving this requires showing 
only upper semicontinuity: F does not exceed F[p] +e in some neighborhood of p. The proof 
involves showing that wavefunctions can be deformed to produce nearby densities without 
costing too much energy. Two aspects of the coarse-grained situation make this possible. 
First, imperfections are hidden. We have to get only the cell averages right. Secondly, X is, 
so to speak, locally finite dimensional. Thus, to show that any required small modification 
can be made in a bounded region £ only requires showing that a finite number of directions 
of modification can be handled. If we had an infinite number of dimensions, as in the 
fine-grained theory, the margin for change in successive directions can shrink to zero. In 
fact, there are infinitely many dimensions — outside E. These are handled in a completely 
different way: the tail of the state can be chopped off, at arbitrarily small energy cost far 
enough out. Then, using the fact that we have mixed states to work with, a new tail can be 
grafted on with the energy cost bounded by Eq. (| 14|) . (For arbitrarily normalized 7, that 
bound extends to (T + V ee ) 1 < -FmaxllpHi, where 7 n- p.) This is the only theorem which 
requires getting behind the densities and working with states, but it is a key ingredient in 
Theorem [42] 

Theorem 13.11 showed that any p £ X^ + is V-representable. Naively, we expect — v[p] to 
coincide with the functional derivative SF/Sp. If — v[p] is in X* (bounded), it is a subgradient 
of F at p. But as discussed in Section [5] it is unclear to what extent directional derivatives 
F[p;6p] agree with it even in that case. Certainly, it is impossible to put all directional 
derivatives together into a linear functional, so that classical Gateaux differentiability is 
out of the question, even if v[p] £ X* . Situations (b) and (c) of Fig. [T] still arise for 
directions Sp which cause immediate exit from X^. But, since those are the only directions 
which cause that problem, it is from the beginning less severe than for the fine-grained 
interpretation. Fortunately, directions Sp which lead immediately out of X^ do not seem 
to hold any physical interest. It is certainly satisfactory if (— v[p],Sp) = F'[p;Sp] for Sp 
satisfying p + sSp £ X£ for some s > 0. That this set of directions really is coextensive 
with 

domF'[p; • ] := {Sp £ X : F'[p; Sp] < +00}, 

follows from 

F [p + sSp] > F[p]+sF'[p;Sp], 

which is a consequence of convexity. It is not ruled out that F'[p;Sp] = —00 for some 
Sp £ dom F'[p; ■]. However, this could happen only if p + s5p falls outside X^ for all s < 0, 
because F'[p; —Sp] > —F'[p; Sp], which is another simple consequence of convexity. 

Theorem 4.2. Suppose p £ X ++ is represented by the potential v. For Sp £ X , either 

a. there is no s > for which p + sSp £ X + , in which case F'[p; Sp] = +oo, or 

b. —00 < F'[p\ Sp] < +00 and F'\p; Sp] = — J vSpdx. 

To paraphrase, if p 6 X^ + is represented by the potential v, then for any Sp £ X we have 
the following dichotomy: Either p + sSp X^ for any s > 0, so that Sp is simply not in 
domF'[p; ■], and F'[p;Sp] = +00, or F'[p;Sp] = — J vSpdr. 

There are domain issues which should be discussed before sketching the method of proof. 
In the fine-grained case, only potentials in X* could be handled systematically, and these 
are bounded by nature: \J—vSpdx\ < \\v\\x» \\Sp\\x- Now, v[p] lies in V, and it is not 
immediately obvious when J —v[p] Spdx is well-defined, since both the positive and negative 
parts can be infinite for some Sp £ X. But in fact, J v[p] ■ dx is unambiguous on all of 
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domF'[p; ■], as the following argument shows. Certainly, v[p]p is integrable. Shift v[p] 
by a (finite) constant to make it positive and split bp into positive and negative parts as 
Sp = (Sp) + — (6p)~, so that —v[p]Sp — — v[p](Sp) + + v[p](Sp)~ . If Sp e domF'[p- ■], then 
{Sp)~ lp is bounded, and since v[p]p is integrable, v[p](Sp)~ is also. Thus, for such Sp, 
J —v[p] Spdr is either real, or — oo. The second possibility is not particularly exotic if v[p] 
diverges as r — > oo. 

The argument used to prove the theorem is similar to the one suggested in Section [2] 
for Gateaux differentiability of F at V-representable densities and which was shown to 
fail in that context. It follows from convexity of F that F'[p\-\ satisfies F'[p;x + y] < 
F'[p;x] + F'[p;y], and F'[p;Xx] — XF'[p;x] for A > 0. That is, it is sublinear, but not 
necessarily continuous. Since J v[p] ■ dx is finite and linear on density perturbations with 
bounded support (call this set JC), the sum functional F'[p\ ■]+ J v\p\ ■ dx is also sublinear 
on such perturbations. Now, the variational principle guarantees that this functional is 
non- negative, so the problem is to show that it is exactly zero on dom F'[p; ■]. If there 
is a direction in K, along which it is not, then on that one-dimensional space, there is a 
nonzero linear functional dominated by F'[p;-] + J v[p] ■ dx. One version of the Hahn- 
Banach theorem says that such a linear functional can be extended to a linear functional 
w on the entire vector space K, which is still dominated by F'[p; ■] + J v[p] ■ dx, since the 
latter is sublinear. In the coarse-grained setting, w is guaranteed to have the same form as 
a potential (in the fine-grained case, distributions, among other things, might arise as linear 
functionals) , and that would make v + w another potential having p as ground state density, 
violating the Hohenbcrg-Kohn theorem. The result is then extended to all of dom F'[p;-] 
using Theorem 14.11 and convexity of F. 

Since all densities are V-representable in the coarse-grained model, a natural next question 
is whether v[p] is a continous function of p. This would say that, if \\p' — p\\i is small 
enough, v[p'] is 'close' to v[p), which would seem to require a topology on V to make 'close' 
meaningful. One topology to consider is the product topology, in which a neighborhood of v 
consists of all v' <E V which are close to v on a specified bounded region, but unconstrained 
outside it, so that open sets are unions of sets of the form 

U(v, E, e)={i)'eV: \v'(x) - v(x)\ < e, \Jx e E}, 

for bounded E. This topology really is weak in the current context: for example, if v n is zero 
for |aj| < n, but goes below — n and above n somewhere outside that radius, it converges to 
zero in the product topology. Be that as it may, we can prove that p i— > v[p] is continuous 
from Xpj with the L 1 topology (as always) to V with the product topology. What needs 
to be shown is that the restriction v[p] fs of v[p] to a bounded region E, viewed simply as 
a vector in a finite-dimensional Euclidean space, is continuous as a function of p. The key 
is to view F[p] as a family F [p p Ie ] of functions of the finite dimensional variable pi's 
parametrized by pts<=, where the superscript 'c' indicates a complement, i.e., E c = IR 3 \ E. 
These are differentiable convex functions on a finite-dimensional space, continuous with 
respect to the parameter p fs , and the very nice properties of convex functions on finite- 
dimensional spaces imply that the derivatives, which are v[p] \s, are continuous in both pfs 
and p\s<=. 

A stronger statement is available, but not, as might be expected, by using a stronger 
topology on V. The crucial observation is that, although v[p] may be unbounded, the 
product v[p]p is always relatively tame. It is integrable, that is, it has finite L 1 norm. 

Theorem 4.3. The map p t— > v[p]p of X ++ into X is continuous with respect to L x -norm. 

Theorem 14.31 savs that the map p H> v[p]p of X^ + into X is continuous with respect to 
L 1 norm. Technically, this is different than saying that p i— ► v[p] is continuous, but it has a 
similar import. The theorem represents a strengthening of the earlier result because, for a 
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bounded region £, 

/ \(v[p']-v[p])\p'dx< f \v[p']p'-v[p]p\ dx + f \v[p](p'-p)\ dx. 

For p' in a small enough neighborhood of p (depending on p indirectly through v[p\) the 
second integral on the right-hand side can be made small and p' will be bounded uni- 
formly away from zero on £, so that some multiple of J s \v\p'\ — v[p] \ dx is bounded by 
/ s — v[p]p\ dx plus a small correction. 

The Proposition in Appendix [Cl shows that this quasi-continuity of v[p] does not extend 
to even the JT*-representable fine-grained densities in Jn ■ 



B. deferred proofs 

The next lemma is preparation for proving Thm. 14.11 

Lemma 4.1. Suppose 70 1— > po is a ground state of T + V ee + v, and let £ be a bounded 
region consisting of entire cells. Then, there exists a neighborhood o/pots; an d a family 
of states 7[pI~e] indexed by p\^€ O such that: 7[pofs] = 7o? the density o/7[pfs] is p\t, on 
£ and po fs c outside it, and Tr{(T + V ee )j[p \^]} is an infinitely differ entiable function of 
Pis- 

Remark. Since F[p] is lower semicontinuous, the lemma implies that F[p] varies continuously 
with density perturbations in a spatially bounded region. It is a key step toward Theorem 
14. li which allows perturbations with unbounded support. 

Proof. First, assume that 7 is a pure state The general case will follow very easily from 
that. 

Let Ur denote the region of iV-particle configuration space where all particles are in the 
interior of a cell R in S. Ur is open. 

Now, we appeal to the unique continuation principle which assures us that \& is not 
identically zero on Ur. This is where we need \& to be an eigenstate of T + V ee + v. The 
simplest version of the principle, only requiring a locally bounded potential suffices here. See 
the Appendix to §XIII.13 in Ref. [3y|. Choose a smooth (C°°) non-negative function gR(x) 
on configuration space, compactly supported in Ur, normalized so that J gR(x)\^\ 2 dx = 1, 
and define ip s (R) = (1 + s{R)gR) 1 / 2 ^! , for a real parameter s(R). Then, 1 + s(R)gR > 
for s(R) sufficiently small, say < e^. Using the Cauchy-Schwarz inequality, the 

finiteness of the various pieces of the energy of 'J, and the fact that derivatives of gR are 
continuous and compactly supported, it is easy to see that || Vi/^ij) || 2 and (4 , s(R)\^ee\' l Ps(R)} 
are smooth functions of s(R). The (coarse-grained!) density of V's(ii) is identical to that of 
^ except in cell R, where it is p(r) + s(R). Now, to complete the construction, just repeat 
with the other cells in S, to get ip s = JJ Re ^{l + s(R)gR) 1 ^ 2 ty. We observe in passing that 
to get smoothness for variations in arbitrary directions requires bounds on the derivatives 
uniform with respect to cell indices. This is a major part of the reason £ must be finite. 

If 7 is a mixed state, we can perform the modification on just one of the pure states in 
its canonical decomposition. □ 

Proof of Theorem 14.11 As F is lower semicontinuous, we only need to show upper semi- 
continuity at p. Further, density can always be added with a density matrix corresponding to 
the desired extra density at an intrinsic energy cost bounded according to F[p] < ||p|| i-F m ax, 
the appropriate variant of Eq. (fT4|) . 

So, we only need to show that, given e, there is 5 > such that F[p'\ < F[p] + e whenever 
p' is in the 5-ball centered at p, Bs(p) : = {p' € X + : \\p' — p\\i < 6} and satisfies p' < p 
everywhere. We will refer to Bs(p) n {p 1 : p' < p} as the 'lower half of Bg(p). 
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Suppose 7 = 2 ■VGV'aXV'al is a ground state with density p. We modify it as follows. Let 
f{x) be a smooth, monotonically decreasing, function M + — s- M + which is 1 for < x < 1 
and for x > 2, and define 3>L: = IIi=i yG 2 ^!/-^)- The modified state is 

Then jl ^ Pl, where pl — outside a sphere S2L of radius 2L. 

As L — > 00, pi converges to p uniformly on any given bounded region £ and the intrinsic 
energy of tends to that of 7. Find L large enough that 

||p-Pi||i <e/(3F max ), 

and 

F[ PL ] < (T + V ee ) lL <F[p] + e/3. 

By Lemma |4.1[ there is some S such that, for p in Bs(pl) and supported in S2L, 

F[p] < F[p L ] + e/3 < F[p] + 2e/3. 

Now, take p' in the lower half of B$(p), and decompose it as p' = mm(p' , pl) + p", so 
that p" < p — pl- Then, according to the previous paragraph, 

F[mm(p',p L )} <F[p] + 2e/3, 

and p" can be added with intrinsic energy cost not exceeding e/3. Thus, F[p'] < F[p] + e. 
□ 

The next lemma is preparation for Thm. !4.2l on directional derivatives of F, and deals with 
the special case of perturbations with spatially bounded support. The subset of X consisting 
of simple junctions, functions which are nonzero only on a bounded set, is denoted /C. 

Lemma 4.2. If p G X ++ is represented by the potential v and bp € K,, then 

F'{p; Sp] — — J vSpdr. 

Proof. It follows from convexity of F that F'[p; ■ ] is a sublinear functional on K,, i.e., F'[p; x+ 
y] < F'[P: x] +F'[p; y], and F'[p; Xx] — XF'[p; x] for A > 0. Since (v, • ) is a linear functional, 
the sum F'[p; •] + («,•) is also sublinear. Since p is the ground-state density of v, F'[p; ■ ] + 
(v, •) > 0. We are trying to show this is an equality. Suppose not. Then, there is some 
Sp € /C such that 

F'[p;Sp] + (v,Sp)>0. 

On the one-dimensional subspace spanned by Sp, there is therefore a non-zero linear func- 
tional A such that X < F'[p; ■] + (v, ■). A version of the Hahn-Banach theorem now implies 
the existence of an extension of A, which we continue to denote by A, such that 

A(.)<F'[p;.] + vV) 

on all of IC. 

The immediate object is to show that A acting on Sp £ JC is represented as 



X(Sp) = / —w(x) Sp(x) dx. 
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FIG. 2. Schematic of construction in proof of Theorem 14.21 



This is not automatic. We have only shown that A is a linear fuctional, not that it is 
continuous. But, for density perturbations nonzero only on a bounded set S, consisting 
of a finite number of cells of there must be such a w, since then the functional is on a 
finite-dimensional space. Such w's for different E's must agree on overlaps and they can be 
patched together to yield a single ^-measurable function. So, the w representation holds 
on all of JC. 
Thus, 



< F'[p;-] + (v + w,-). 

This relation appears to say that p is the ground-state density for v + w. If w is bounded 
below, that is correct and we can therefore apply the coarse-grained Hohenberg-Kohn the- 
orem to conclude that w = 0. But if w is not bounded below, that argument is not 
immediately applicable. In that case, consider, the density perturbation SpR which just 
adds 1 to some arbitrary cell R. Since F is bounded above and below, and is convex, 
— w(R) — v (R) < F'[p; Sp R ] < F[p + Spn] — F[p] < F max . Therefore, since R is arbitrary, 
v + w is bounded below, so the Hohenberg-Kohn argument actually does apply to v + w, 
after all, showing that w = 0. □ 

Proof of Theorem l4.2l Case (a) is clear. So suppose there is s > such that p+s5p £ X + . 
Renormalizing Sp if necessary, assume without loss that p + sSp G X + for s < 1. 
Fix e > 0. Using Theorem 14. 11 find S such that 



F[p']-F[p] <e/4, p'eB s (p), 



and then a bounded region S large enough that ||<5p[£ c ||i < 5. 
Now, we split F[p + s 5p] — F[p] as 



(16) 



F[p + s Sp] - F[p] = [F[p + s Sp] -F[p + s 6p\v] 

+ (p[p + sSph}-F[p]), (17) 

and deal with the two terms ("1st line" and "2nd line") on the right-hand side following a 
strategy illustrated by Figure [5] 

For the 2nd line, since Sp\s has bounded support, Lemma 14.21 savs that 

F[p + s6p\z] -F[p] = -s / v5p\z dr + o(s). 



As for the 1st line, convexity of F implies that 

F[p + sSp] < (l-s)F[p + s5p\ s ] + sF[p + s5p\ s +5p\v>], 
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so that the 1st line satisfies 

F\p + a8 P ]-F\p + s6p\v]<a(F\p + 88p\v +6p\&] - F[p + s Spfc]) ■ (18) 

For s small enough (Fig. EJ, both p + s5p\s and p + sSp\s +Sp\^c are in Bs(p), so that 
their intrinsic energies do not differ by more than e/2. Injecting that fact into inequality 

F[p + s5p}- F[p + s 5p\x] < se/2. 

Putting everything back into Eq. (fTTl) . F[p + sSp] — F[p] < s (e/2 — J v5p\z dx) + o(s), 
showing that 

F'[p;5p\ <e/2- / v5p\ s dx. 

In the limit E /• the integral tends to J v 8p dr whether or not 8p is in the domain of v. 
If it is not, the value of the integral is unambiguously — oo. Taking now the limit e — > gives 
F'[p;5p] < — JvSpdx. The variational principle already secured the opposite inequality, 
thus we have equality: F'[p; 5p] — — J v 5pdx. □ 

The final lemma-theorem pair in this section is aimed at getting something resembling 
continuity of v[p] as a function of p. 

Lemma 4.3 ('Local continuity' of potential). The potential v[p] is a continuous function 
of p with respect to the L 1 topology on X ++ and the product topology on V. 

Proof. To say that v[p] is continuous with respect to the product topology is to say that 
v[p] \s is a continuous function of p. Thus, fix a bounded region S and view 

F[p] =F[p\x+ph4 

as a function of the two variables pi's and p\sc. 

For fixed p\z e , F[p] is a bounded convex differentiable function (according to Lemma l4.2p 
of p\-£. Temporarily denote the derivative with respect to p\y, by 

D s F\p]. 

According to a theorem of convex analysis (Appendix IB1 (6)), boundedness of F and finite- 
dimensionality of the variable p\s imply that D^F[p] is a continuous function of p\^. 

On the other hand, for fixed p\%, F[p] is a continuous function of p according to 
Theorem 14.11 Again from boundedness and finite-dimensionality conclude that D^F[p] is 
also a continuous function of p fs=. For, if not, there would be a sequence of values r„ 
converging to p\^a such that + r n ] converges to F[p\^ + p\sc], but D^F[p\^ + r„] 

does not converge to D^F[p\^, + p\s c ]- According to Appendix [B] (5), that cannot happen. 

But, the derivative of F[p] with respect to p \s is precisely v[p] fs- So the previous 
statement says that v [p]\s is continuous with respect to both pfs and pts^, which is to say, 
it is continuous with respect to p. □ 

Proof of Theorem 14.31 Suppose not. Then, there exists e > and a sequence p m — > p 
such that 

\\v[p m }p m - v\p]p\\i > e 5 for all m. (19) 
We will derive a contradiction. 
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According to our convention [Eq. (fT5|) ] on the constant offset of v[p], 

/ v{p m ]p m dx = -F[p m ] - / v[p m ]p m dx, (20) 

for any bounded region E. As m — > oo, the first term on the right-hand side converges to 
F[p] by Theorem 14. II Also, by Lemma FOl the integral on the right-hand side converges to 
J s pete, since the integration is over a finite region. Actually, that Lemma shows more. 
The contribution to Inequality (|19p from integration over E tends to zero as m — > oo. 

So, as to — y oo , both terms on the right-hand side of Eq. ([20]) tend to their counterparts 
with p m replaced by p. Since the entire equation also holds with that substitution, the 
left-hand side must converge to L c v[p] pdx. Abbreviating !„(£) = J EC v+/? m dx, /„(£) = 
J sc v^pra dx, and similarly with no m subscript, that means 

J+(E) -/-(£)-> !+(£)- I" (E). 

Now, both I + (E) and /"(E) tend to zero as S / M 3 . This shows that 

J+(E) -> I m (E), as S / M 3 and m -> oo. 

On the other hand, as has been remarked, inequality (fT!?]) is almost entirely carried by E c 
as to — >• oo, for any E. So, for large enough E and to, 

/+(£)+ 1- (£)>e/2. 

The last two displays together show that, for large enough E and to, both i^(£) and 
i^(£) are bounded away from zero. But that is clearly impossible. v m is bounded below 
(v~ < F max ), so that 

p m dx. 



Yet, Jy, c p m da: tends to zero as E /* M 3 and to — > oo, since p m — >p. This contradiction 
finishes the proof. □ 



5. KOHN-SHAM THEORY 



This section concerns the implementation of Kohn-Sham theory within the coarse-grained 
model. Since all densities in X~^ + are both interacting and non-interacting V-representable, 
the endeavor is off to a good start. Two additional issues are a little subtlety in the definition 
of the Hartree energy, and the relation v xc = SE xc /5p. 

Of course, the Kohn-Sham decomposition of the intrinsic energy for p G X + looks just 
like that for the fine-grained theory, Eq. JTJ]): F[p] = T s [p] + E H [p] + E xc [p\. Only the 
interpretation is changed. E xc [p\ is defined only on X^, since outside, both F[p] and T s [p] 
are +00. Similarly, Eq. (|13p is taken over and used to define the exchange-correlation 
potential, as v xc [p] := v s [p] —v[p] — <j>[p] on X^ + ■ 

The definition of the Hartree energy requires supplementation, because a coarse-grained 
density only fixes the total particle number in each cell. A choice must be made here which 
does not arise in the fine-grained theory, and there seem to be two possibilities. Different 
choices here would imply slightly different E xc . One possibility is to use a density which 
is uniform throughout each cell. In that case, we would begin using this surrogate for a 
coarse-grained equivalence class in a more explicit manner, and the Hartee energy would 
depend only on the coarse-grained equivalence class. The second possibility is to use one of 
the fine-grained densities which minimize the intrinsic energy. This certainly seems natural 
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in some ways, but it raises problems. First, which intrinsic energy, F or T s ? Even if that is 
decided, say for F, there may be multiple fine-grained densities which minimize F, calling 
for another choice. One of the nice features of the Hartree energy is its explicitness. That 
would be lost with this choice, and possibly the ability to establish continuity of Eh as well. 
For all these reasons, we will take the first choice: the Hartree energy is calculated according 
to a uniform distribution of charge in cells. 

Off X + , p can be negative, and if we interpret that as positive electrical charge, Ejj 
becomes well-defined on all of X. It is also[37] convex, bounded below by zero, and bounded 
above by a multiple of \\p\\i, the multiple being determined by the energy to charge a single 
cell. These properties imply that Eh is continuous as a function of p G X. 

The derivative of Eh, which will be denoted <f> is explicitly computable as an element of 
X* . From ip(x) = e 2 J p(y)/\x — y\ dy, <fi is obtained by averaging over cells. The map 
p i y cf>[p] is linear, and it is not difficult to see that it is continuous from X to X* , which 
is to say, \\<j>\\oo < C||/°l|i- Thus, is a Frechet derivative of Eh- The Hartree potential 
behaves better than does v[p\. 

Now we turn to the exchange-correlation potential. It follows immediately from the 
definition that p > v xc [p}p is ./^-continuous, since the other three potentials have this 
property. Is v xc the derivative of E xc ? More precisely, does it coincide with the directional 
derivative E'[p;-}7 Consider what can go wrong. If Sp is such that p + s5p is not in XZ 
for any s > 0, then E' xc [p; Sp] has no value, finite or infinite, because E xc is not defined off 
X^. Otherwise, we should have E' xc [p; Sp] — F'[p; Sp] — T' s [p; Sp] — E' H [p; Sp] = (v[p],Sp) — 
(v s [p],Sp) — {<j)[p],8p). This is correct as long as at most one term on the right hand side 
is infinite. Either F'[p;Sp] or T s [p;Sp] might be — oo. A condition which can be imposed 
directly on p to make sure E' xc [p; Sp] and (v xc , Sp) both exist and are equal is that p + s Sp is 
in Xjf for s in some open interval around zero. The fairly solid status of the coarse-grained 
exchange-correlation potential is in stark contrast to its fine-grained counterpart. 

6. MULTISCALE MODEL AND LIMITS OF ZERO COARSE-GRAINING SCALE 

The previous section aimed to show that coarse-graining cures some of the bad behaviors 
of the continuum theory. But not everything about the continuum theory is bad, and it is 
desirable that coarse-grained models be good approximations to it in certain respects. This 
immediately leads us to ask about continuum limits. What happens as the coarse-graining 
scale is taken to zero? Is the continuum theory recoverd as smoothly as possible? Are there 
misleading limits? These are the concerns of the present section. 

A. ideas 

A straightforward way to approach a continuum limit is to use a sequence of ever-finer ac- 
ceptable partitions n = 1, 2, 3, ... of the sort introduced in Section [31 with the maximum 
cell diameter of denoted D n , tending to zero as n — > oo. For the sake of the Poincare 
inequality in Thm. 16.11 we impose the technical condition that the cells be convex. Finally, 
in order to have the collections of coarse-grained densities strictly increasing with level, we 
require the cells of Vpn+i to be obtained by subpartitioning those of For example, these 
requirements are satisfied if we take the cells of to be those of a simple triclinic lattice 
generated by lattice vectors 2~ n dk with a%, 02, and 03 noncolinear. 

As in Section^ each gives rise to a space X n of equivalence classes of densities. The 
object now is to make contact with the theory sketched in Section [21 where X — L 1 n L 3 
and its subset Jm figured prominently, so we will equip X n with the L 1 n L? norm as well, 
instead of the L 1 norm as we did with the single-scale model. The point is that this does 
not require any essential change to what was done in the previous two sections because for 
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any fixed X n , the L 1 n L 3 and L 1 norms are equivalent: ||/||3 < c n ||/||i for / G Af n . For 
a sequence on coarse-graining scales tending to zero, though, convergence with respect to 
L 1 n L 3 is more stringent than L 1 convergence. In studying the single-scale model, it was 
convenient and harmless to represent elements of X by ^-measurable functions (that is, 
constant on cells). In the multiscale setting, this conflation is not so innocuous. Therefore, 
we introduce the injective isometry 

i n : X n ^X, (21) 

which takes p £ X n to the unique ^n-measurable function belonging to the equivalence 
class in X associated to p. 

For any n, each g € X belongs to one equivalence class corresponding to an element of 
X n , which will be denoted by ir n g. As a notational cue, elements of X will generally be 
denoted by a rho with tucked- in tail (g), and elements of one of the X n by a normal rho 
(p). We refer to ir n g as the (scale-n) projection of g. This nomenclature is doubly justified, 
since i n X n is a closed subspace of X and L n ir n g is precisely the L 2 orthogonal projection 
of g onto this subspace. We wish to think of coarse-grained densities as formal objects 
independent of X , and thus stop short of simply identifying them with equivalence classes 
in X . Consistently with the notation just introduced, the equivalence class corresponding 

to p is 7T~ /?. 

Each X n has its own intrinsic energy functional F n . For p € X n , 

F n [p] :=M{F[q] -.gen-ip}. 

As a notational convenience, we write F n [g] for F n [ir n g). As discussed in Section [3j if 
p € (X n )~jy, the infimum is attained at some fine-grained density, and possibly more than 
one, though such degeneracy is not generally expected. For p € (X n )j^, the entire set of 
such minimizers will be denoted 

A n p:={Q€Tr- 1 p:F[ Q }=F n [p}}. (22) 

In particular, if p is misnormalized or somewhere negative so that -F™[p] = +00, then 
A n p : = 0. Thus, A n is a set-valued function or multifunction\27\. l38l. l39j . 
The collection of all coarse-grained densities at all scales is denoted 

00 

X°° := (J X n , (23) 

n=l 

and for p e X°° , the scale function is defined by p G ^ scale (p). The following chain of 
inclusions, all of them strict, then holds: 

l 1 X 1 Cl 2 X 2 C ••■ C lX°° =X\ J N CX. (24) 

On p € X°° , 1 acts as t sca | e (p). The scale index on t (and on A) is strictly unnecessary, but 
will be written at least when disambiguation is thereby provided. Note that i n X n is a closed 
vector subspace of X, but iX°° is not a vector space at all. It does not make sense to take a 
linear combination of coarse-grained densities at different scales. The fine-grained densities 
of interest are in Jn . 

Turning to potentials, we have a chain of spaces 

00 

V 1 CV 2 C •■■ C V°° := (J V", (25) 

n=l 

where V" consists of *P„-measurable functions bounded below. The coarse-grained theory 
of previous sections provides a representing potential map v n : (X n )~^ + — > V™. From the 
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perspective of the fine-grained theory, v n [p] is the representing potential of A n p, and in 
making contact with that theory, the subspace (X n )* — {v e V™ : < 00} of the 

potentials with finite L°° PI L?l 2 norm, and understood to be carrying that norm, is of more 
interest. These nest as 

(X 1 )* d(X 2 )* C ■ ■ • C X*°° := U~ =1 (A- n )* C X* . (26) 

Note that for the potentials, there are no embedding maps analogous to 1. A coarse-grained 
potential really is a *P„-measurable function. 

Now we can begin to investigate how these ideas fit together. The first set of questions 
involve convergence, or lack thereof, of LTr n g and Air n g to g. The positive answers hold 
implications for convergence of F n [n n g] . 

Theorem 6.1. Given g in Jn , i^nQ Q- F n [o\ *s increasing with n, and F n [g] = F[A n g] 
F[g] as n — >• 00. Both Air n g and iTr n g converge to g in L p norm as n — > 00 for 1 < p < 3, 

with \\g - ATT n g\\ p ,\\g - in n g\\ p < cD^ p)/2 F[g}( 1+ Py 4 . 

Convergence ||wr n £ — g\\x — > is just a general property of LP functions under cell- 
averaging. L 2 convergence is immediate from the description of 7r n as an orthogonal pro- 
jection. Convergence in LP for other p is proven by Jensen's inequality. Although it is 
not exactly the sort of approximation coarse-graining was designed to produce, it is not 
surprising that Ln n g tends to g in X as n — > 00. It is much less obvious that Lir n g should 
converge to A-K n g since Air n g is a moving target. We might just as well ask how spread out 
7r J 7 1 p is, as a subset of X. Every element g e n~ 1 p must match ip down to scale scale(p), 
but might differ arbitrarily much at "subgrid" scales. However, there are limits imposed 
by F[g] through the lower bound in ([5]). Roughness at subgrid scales come at the cost of 
intrinsic energy, as intuition suggests. Using the bound with a Poincare inequality enables 
us to prove a bound on the L p norm \\iTr n g — g\\, p for 1 < p < 3 in terms of F[g]: 

||wr we - Q\\ p < c p Dl(F[g\/Dl)^'\ 1 < p < 3, (27) 

where D n is the minimum cell diameter in Thus, the set of fine-grained densities in 
ir~ 1 p with low intrinsic energy, ir~ 1 p n {F < M} is bounded in LP norm. Although the 
bound in (|2"T|) holds for p = 3, it does not show convergence in that case since the factors of 
D n disappear. This clears the way to show that hir n g — > g in L p . For, it is a triviality that 
F 1 ^} ^ F 2 \o\ < F 3 \q\ < ' ' ' < F[q]> since each successive term represents a minimization 
with additional constraints (over a smaller set). Since F[Air n g] = F n [g], there is a common 
bound for both \\g — Lir n g\\ p and ||A7r„p — t7r„^|| p improving with n. We have noted that 
is increasing with n toward F[g], but could conceivably fail to reach F[g] in the limit. 
But, since F is L 1 lower semicontinuous and A n g converges to g in L , the gap does close. 
No intrisic energy somehow goes unaccounted for in the limit. 

We can broaden the scope of the above a little by considering not only sequences going 
"straight up" the hierarchy of spaces along 7r„g, but also those converging toward that one. 
A sequence (pj)f^i in X°° is said to converge to g G Jm ii \\g — Lpj\\x — > as j —> 00. Since 
iir n g is the L 2 projection of g onto tX n and is at nonzero L 2 distance from it, g is at nonzero 
A-distance from iX n , Thus, norm convergence of ipj to g implies that scale(pj) — > 00. And 
therefore, since Thm. 16. II has shown that iTT n g — > g, the condition defining "coarse-grained 
sequence converging to g" can equivalently be written as \\i(pj ~ ^ ni g)\\x — > 0, where 
rij = scale(pj). 

The convergence results of Thm. I6.1l can thus be generalized mildly to sequences (pj)j°. 1 C 
X°° converging to g, but the bounds on F nj [pj] (rij := scale(/5j)) are not automatic in that 
case and must be imposed as hypotheses. 

Corollary 6.1. Let (pj)^_ 1 be a coarse-grained sequence converging to g. Lf Fj := i^Apj] < 
M < 00, then Apj — »• g in L p , where rij = scale(pj). 7/limsupFj < F[g], then Fj — ► F[g\. 
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We now shift attention to convergence questions related to representing potentials. 

Theorem 6.2. g G Jn is X* -representable if and only if there exists a coarse-grained 
sequence (pj)f%i converging to g such that {\\v nj [pj] \\x* ■ 1 < j < oo} is bounded (rij : — 
scale(/3j)J. In that case, v nj [pj] — > v[g\ in weak-* sense and F nj [pj] — > F[g\. Further, there 
exists such a sequence such that the convergence v nj [pj] — » v[g] is in X* -norm. 

Unfortunately, the tools do not seem to be at hand to say anything significant beyond 
the context of A*-representability. However, in the latter context the closedness of the 
graph of dF and the Ekeland variational principle, both discussed in Section^ are powerful 
tools. Suppose (pj)j? =1 is a coarse-grained sequence converging to g € JTn, with representing 
potentials v 3 = v[Apj}. If (v 3 )j°. 1 is merely bounded in A*-norm, then g is A*-representable 
and ir* — > v[g] in the weak-* topology. At first, this might look like a straight transcription 
of the closedness of graph dF discussed in Section[2l If it were the case that A nj pj converged 
to g in X, that would be correct. For, bounded subsets of X* are weak-* compact, and any 
weak-* cluster point v would be in dF[g\. Since the Hohenberg-Kohn theorem guarantees 
that d-F 1 ^] is a singleton modulo constants, the entire sequence would have to converge. 
However, it does not follow from ipj — > g that A nj pj — > g. This circumstance requires using 
the Ekeland variational principle. Showing that g is almost a ground state of v 3 means that 
the sequence (A nj pj , v 3 )"^° can be replaced by a new one (pj^ 3 )^ such that ipj — >■ g and 
\\v 3 — v 3 \\x* — > 0. Closedness of graph OF then applies in the ordinary way. 

In the other direction, suppose that g G Jn is A*-representable. Then, according to the 
theorem, this is accurately reflected in the coarse-grained hierarchy: there is a coarse-grained 
sequence (pj)°^ 1 converging to g, such that v 3 — > v[g] (in this case we get norm convergence 
of the potentials). The proof of this direction is perhaps a little more interesting. To prove 
this, we find an approximation to v[g] in (X n )* for which g is nearly a ground state. Then 
we apply the Ekeland variational principle, but in X not in X . This produces a sequence 
of pairs (pj,v 3 )^° with the desired property. The disappointing aspect of this half of the 
theorem is that we cannot show that v 3 [ir nj g] — > v[g\. But, there was no good reason to 
suppose that true, anyway. 

We finish this section with some reflection related to the computational and physical 
appropriateness and significance of topologies on X°°. In preparation, we have a grab-bag 
theorem about the multifunction A„ . A remarkable aspect of the bound (|2~T|) is that it can 
be used to show that A„ preserves compactness. 

Theorem 6.3. A n takes values in convex, L 1 -compact sets. If K is ^-compact in X n , 
then A n K is L 1 -compact. A n has a (L x L 1 ) closed graph and is upper semicontinuous 
(meaning, given p and e > 0, there is there is 8 > such that — p\\i < 5 implies that all 
of A n p' is within e of A n p.) 

Proof. That A n p is convex is trivial. That it is L 1 compact is a special case of the compact- 
ness of A n K, shown below. Closedness of the graph of A„ is a simple consequence of the 
L 1 lower semicontinuity of F. For, suppose {pj,gj)j^i is an L 1 x L 1 Cauchy sequence in 
graphA„. Then Pj ^ p e {X n )+ and gj g € X+. Also, F[g] < limj^ F[gj} = F n [p), 
where the inequality follows from lower semicontinuity of F and the equality from continu- 
ity of F n . Since g € tt^P, F[g] > F n [p], so that actually F[g] — F n [p]. In other words, 
g G A n p, as was to be shown. 

To see that A n preserves compactness, let K C (X n )~^ be L 1 compact. Since F n is 
continuous, it is bounded above on K, say by c. Lemma 16.31 in the next subsection and the 
previous paragraph then show A n K to be a closed totally bounded set, hence compact. 

Suppose upper semicontinuity failed at p. Then, there would be a sequence pj — > p and 
gj G A n p such that every gj was at L 1 distance greater than e from A n p. But, {pj : 1 < j < 
oo}U{p} is compact. So {gj : 1 < j < oo} is relatively compact, and therefore by closedness 
of A n contains a subsequence converging to an element of A n p. The contradiction proves 
upper semicontinuity. □ 
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This theorem shows that, for compact sets K C X n , not only is A n K not very far from 
K (Theorem 16. ip . it is also not much larger. 

There are a couple of ways to look at the a single space X n , the L 1 or L 1 n L 3 metrics 
are pullbacks via i from X. 

One way to look at the L 1 or L 1 n L 3 metrics on a single space X n is to recognize 
that they correspond to the same topology, namely the weakest topology which makes the 
particle number in each cell of and the total particle number continuous. This topology 
seems natural when considering coarse-grained densities. Alternatively, we recognize them 
as pullbacks via u or tt^ 1 from X. While the connection with i is fairly obvious, it is also 
the case that 

A (P)P') : = IIPjP'IIi = mf {ll £ - e'lli : Q € \ € 7r~ 1 p / } ) 

and similarly for the X (L 1 n £ 3 ) norm. That is, the distance is just the distance between 
the sets (equivalence classes) w^p and Tr^p'. 

On the other hand, we have some tendency to see p also as a kind of surrogate for the 
low-intrinsic-energy densities in ir~ 1 p, specifically A n p. After all, that is what the definition 
of F n was all about. This inclines us to consider the metric 

d ( £\p,p') = Lnf{||e-e'||i : ge A n p, g e A„p'} 

on (X n )~p f . (Only (X n )^- can be metrized this way since A is null-valued outside.) It is clear 
that this is at least as strong as the d\ metric. What is interesting is that it is topologically 
equivalent to d^ . In other words, d£ (p, p') is continuous with respect to the d\ topology. 
This conclusion follows from Thm. 16.31 

Moving from the single-scale of X n to the multi-scale setting of X°° , uncovers a new 
oddity. Earlier in this section, for purposes of establishing the coarse-grained model as a 
good approximation of the fine-grained theory, distances between coarse-grained densities 
at different scales were computed as the A-norm distance between their images under l or 
7T _1 . Arguably, this is somewhat too crude. Certainly, it results in some coarse-grained 
densities at different scales being at zero distance from one another, which does not seem a 

(n) 

desirable outcome. Likely, the analogous extension of d A does not have that problem, but 
the uncertainty just points to a different difficulty, which is that of computing d\. 

From a computational perspective, a coarse-grained density is a description with a degree 
of precision related to the scale and it seems that the scale function should be continuous. 
An example of a metric which achieves that is 

d(p,p') = || tp - L.p'\\x + |scale(p) - scale(p')|. (28) 

It should be noted that all the results of Section[3]are true for X°° under this metric. There 
is no mathematical depth whatever to this metrization of X°° , and it is effectively the same 
as just treating the different scales as incomparable. But maybe the mere construction of 
X°° will cure us of the idea that some fundamental length scale is inherent in coarse-graining 
per se. 



B. deferred proofs 

Proof of Theorem 16.11 The convergence of Lir n Q to g in X norm hinges on Jensen's 
inequality. Since ?/ H \y\P is convex for 1 < p < oo, 

W\L edx ) -W\L QPdx ' 
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for each cell il of a partition (volume |0|). Summing the left-hand side over cells yields 
||7rp||P, and summing the right-hand side yields \\g\\^,- Thus, the operator g i— > g — iir n g is L p 
bounded, with bound 2. Since C C (R 3 ) (compactly supported continuous functions) is dense 
in L p , we can find / g C C (R 3 ) with ||/ — g\\ p < e/2, so that || (/ — wr„/) — (g— Ln n g)\\ p < e for 
all n. But, / is uniformly continuous, being compactly supported, so that ||/ — wr n /|| p — > 
as n — > oo. Thus, ||t7r„p — g\\x — > 0. 

A quantitative estimate of || An n g — g\\ p and ||i7r„p — g\\ v hinges on an L 1 Poincarc 40J 
(or Poincare-Wirtinger [4lj ) inequality requiring convexitv|42l l43j . but not regularity, of the 
domain. If / is a function such that it and its (distributional) gradient V/ are integrable 
over the convex bounded region fl with diameter Diam(O) (i.e., / e W 1A (fl)), 



f |/-</) n |daj<£Diam(fi) f \Vf\dx, 
Jo. 1 Jn 



where (/)n is the mean of / over Q. Applying this to g and A7r„p on each cell and summing 
the results, 

J \q- Ln n g\ dx<^D n J \Vg\ dx, (29) 

and similarly for h.ir n g, more precisely (as should be understood for such a locution) for 
each element of K-K n g. To make use of this, we bound ||Vp||i by using inequalities (|5l6p and 
the Cauchy-Schwarz inequality, as 



J \Vg\dx = J 2g 1 / 2 \Vg 1/2 \dx 

(r \ ( r ^ 1//2 

gdxj I / \Vg 1/2 \ 2 dx 



<c'NF[g} 1/2 . (30) 

Therefore, since all the F n [g] = F[Kn n g\ are bounded above by F[g], both \\g — tTr n g\\i and 
\\ATT n g — t.n n g\\i are bounded by cDnFlg] 1 / 2 . We obtain not only convergence, but a bound 
on the rate of convergence. 

To extend this to LP for 1 < p < 3, we make use of a Holder inequality. For /, g g L 1 n L 3 , 

\\f-g\\ P <\\f-g\\t 1)/2 \\f-9\\t p)/2 

<2^ 1 )/ 3 (||/||3 + || 5 ||3)^ 1 )/ 6 || / _. g ||(3-^. 

Since each ||A7r„£>||3 as well as \\g\\3 is bounded by a constant times F[g], combining this 
with L 1 convergence of h.ir n g to g implies L p convergence as well. 

Finally, we prove upward convergence of F n [g] to F[g\. Since represents a mini- 

mization with increasing constraints as n increases, F 1 ^] < F 2 [g] < F 3 [g] < ■■■ < F[g] 
is trivial. On the other hand, as discussed in Section [2l F is L 1 lower semicontinuous. 
Since F n [g] — F[Air n g] and we have shown that Air n g converges to g in L 1 , it follows that 
liminf„F™[£>] > F[g\. Thus, F n [g] /* F[g\. 

□ 

Proof of Theorem ICT2l 

"If" direction: To simplify notation, write nj for scale(/9j) and Vj for v nj [pj\. Also, the 
v nj [pj] given by hypothesis satisfy our convention for fixing the constant offset in represent- 
ing potentials [Eq. (IT5|) ]. but we otherwise lift that convention. It will be restored at the 
end. 

The idea is to show that g nearly attains the ground state energy of Vj for large j, so 
that the Ekeland variational principle can be called in. First, we show that for j large 
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enough, the potential energy of g in Vj is not much more than that of A n . pj . Since Vj is 
^P nj -measurable, it cannot distinguish g from ir nj g or pj from Apj. Thus, 

(vj , g) = (Vj , n n:l g) = (vj , Apj) + (Vj , (tt„ 3 g - Pj )) , 

so that, by hypothesis of boundedness of the potentials, for some M, 

\(vj,g) - (vj,A nj pj)\ < M\\Tr nj g- pj\\x -> 0, as j oo. (31) 

Second, we examine the intrinsic energy. Since A nj pj are the ground state densities of Vj, 
and since as just shown, A nj pj and g look increasingly alike to Vj as j — > oo, F nj [pj] cannot 
be much larger than F[g] for large j, or eventually g would have lower total energy in v J 
than Apj. That is, 

limsupF[A„^] < F[g\. 

Thus, by Corollary EU 

lim F[A njPj ] =F[g}. (32) 

It follows from ([31]) and that, given e, F[g] + (vj,g) < E[vj] + e 2 for j > j(e). 
Appealing to Ekeland's variational principle, conclude that there exists (g e ,v e ) £ graph Si 7, 
satisfying 

\\v e - Vj^Wx* < e and \\g e - g\\x < e. 

Suppose for the moment that Vj converges to v weak-* in X*. In that case, for e — > 
v e — > v weak-*, while g e — > g in X-norm. From the closedness of graph dF, it would follow 
that g is a ground state density of v. 

But, because {vj} is bounded in X* , it is weak-* compact, hence has weak-* accumulation 
points. For any such, say v, we can extract a subsequence of {v j} converging to it weak-* and 
the argument above applies. On the other hand, the Hohenberg-Kohn theorem guarantees 
uniqueness of the representing potential for g. Putting it together, we conclude that the 
original sequence Vj = v nj [pj] converges weak-* to v and g is a ground state density of it. 

Now check that v satisfies the constant-offset convention, Eq. (fTS")) . From weak-* con- 
vergence, (vj,g) — > (v,g). Together with the second displayed equation, this shows that 
(vj, P j) — > (v, g). On the other hand, F n ' [ P j] — > F[g], so (v, g) + F[g] = 0. 

"Only if" direction: Let the tolerance e > be given. Compactly supported continuous 
functions are dense in X* , and any such can be uniformly approximated by some element 
of W^L Q (X n )* , so the latter is also dense in X*. Thus, there exists n and v' e (X n )* with 
\W — v\\x* as small as desired. E[v] is Lipschitz continuous fl3| with respect to AT*-norm, 
so n and v' can be chosen such that F[g] + (v', g) — E[v] + (v' — v, g) is a close to E[v'] as 
desired. In fact, since v 1 £ (X n )* and F n [g) — > F[g] as n —> oo by Theorem 16.11 n — n(e) 
and v' e (X n )* can be chosen so that \\v' — v\\x* < e/2 and 

F n [n n g} + (v',n n g) - E[v'\ < e 2 /A. (33) 

As in the "if" direction, appeal to Ekeland's variational principle again, but this time 
using the pair [X n , (X n )*] together with F n , to conclude to that there exists a pair (p, v) € 
graph dF n C X n x (X n )* such that 

\\v - v\\x* < e 

and 

11/5 - n n g\\x < e. 

Taking a sequence of values of e tending to zero then produces a sequence as described in 
the statement of the Theorem. □ 
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Lemma 6.1. Let K C X n be toally bounded with respect to L 1 D L p for 1 < p < 3. Then 
for any < c < oo ; t^" 1 -^ (~l {F < c} is also totally bounded with respect to L 1 n L p . 

Proof. Since nothing in 7r _1 p has finite intrinsic energy unless p is positive and properly nor- 
malized, we tacitly assume that all coarse-grained densities involved have those properties. 
With e > given, we show the existence of a finite (3e)-net for tt~ x K n {F < c}. According 
to Theorem l6.1[ there is m > n such that \\i> m p~ QW^nL? < e f° r an Y Q S P with F[g] < c. 
Then, every p in Tr^ifn {F < c} is within e of (t of) some point of (tt™)~ 1 K n the 
subset of (<f m ) + which projects to if C A"™. Thus, all that is needed is a finite (2e)-net for 

Now note the following two simple facts. For any scale m, the L p norm of an element of 
{X m ) + is bounded in terms of its L 1 norm. Also, a non-negative element of (-7T™) -1 /? has 
the same L 1 norm over any cell of ^P„ as does p. As a result, since K is totally bounded 
by assumption, there is a bounded *p„-measurable region f2 such that the L 1 n LP norm of 
every element of (tt'^ l )~ 1 K n (A" m ) + over fi c (the complement of Q) is less than e. 

Thus, it will suffice to restrict everything in (ir™)~ 1 K<~) (X rn ) + to VL and find a finite e-net 
for that set. But that is simple because there are a finite number of degrees of freedom in 
those truncated X rn densities and each one is bounded due to total boundedness of K. □ 



7. CONCLUSIONS 

The results of this paper show that coarse-grained models of DFT are mathematically well- 
behaved and are good approximations to the standard, fine-grained, interpretation in that 
non-pathological aspects are faithfully reflected. These properties make the coarse-grained 
models good regularizations of the fine-grained theory. In such a role, coarse-grained models 
provide a controlled arena in which to work, and also shed some light on the pathologies of 
the fine-grained theory. 

Beyond that, however, I have argued that the coarse-grained models provide a superior 
interpretation of DFT. The standard interpretation of DFT takes the somewhat reflex- 
ive view that a density involves an specification of infinitely fine spatial resolution. The 
coarse-grained interpretation understands each density to be specified with limited resolu- 
tion, though there is no limit to that resolution. In DFT, one keeps track of spin density, 
but not degrees of freedom of the state which do not affect it. Those are subject to an 
automatic energetic selection. The step to coarse-grained DFT is very similar, except that 
the degrees of freedom which the formalism gives over to automatic energetic selection are 
the distribution of density within cells — the short length scale degrees of freedom which 
resulted in the fine-grained interpretation not being a good model. Paired in a natural way 
with a reinterpretation of external potential, this yields a model of DFT where everything 
works the way it is supposed to. 
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Appendix A: functional analysis survival kit 

This Appendix contains a quick review of the basic functional analysis used in this paper. 
Readers who need a quick reminder of a definition or notation should find what they need. 
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With luck, readers unfamiliar with functional analysis may find enough to follow the rough 
outlines of the ideas and arguments. Everything here is standard. 

Topology and metric spaces. Questions of continuity and convergence are a major 
preoccupation in this paper. The structure on a set S which makes these notions meaningful 
is topology. Neighborhood bases are a convenient way to specify a topology. A neighborhood 
Af x base for x G S is a collection {U a } of subsets of 5* containing x, and satisfying the 
condition that, whenever U a ,Up G N x , then there is some ?7 7 G N x such that XJ 1 C U a , UpC] 
N x . Any superset of a member of N x is a neighborhood oi x k subset T of S is open if, 
whenever T contains a point y, it contains an entire neighborhood of y. A sequence (i n )J° =1 
in topological space S converges to x if, for any U G N x , the entire tail of the sequence from 
some n (depending on U) onward is contained in U. A subset F G S is said to be closed 
precisely when its complement S \ F is open. The interior of a set T is the largest open set 
contained within T and the closure of T is the smallest closed set containing T. 

When one refers to "the topology T" it means the collection of open sets. One topology 
T' is called stronger than another, T, if T C T' ■ For norm topologies (see Banach spaces 
below), the condition can be expressed by || • || < c|| • ||' for some c > 0. 

Metric spaces are a particularly pleasant sort of topological space. One choice of a neigh- 
borhood base of i in a metric space (S, d) is the collection of open balls B(x,r) — {y G 
S : d(x, y) < r} for rational values of r. Closed balls, or all real values of r would serve 
equally; they all specify the same topology. The metric space S is even nicer if it is complete, 
meaning that if the sequence (x n )^Li is a Cauchy sequence, then it actually converges (to 
some x G S). 

T C S is compact if it has the property that for any covering by open sets, T C U a U a , 
there is a finite subfamily U ai , . . . U an which still covers T. Roughly speaking, a compact 
set is almost topologically finite. If S is a metric space, then T is compact if and only if 
it is complete and totally bounded. The latter condition means that, given e > 0, T can be 
covered by some finite collection of balls of radius e. For the metric space T, compactness 
is equivalent to sequential compactness which means that every sequence has a convergent 
subsequence (to a point in T). 

Banach spaces. Let V be a normed vector space, with norm || • ||. The norm provides 
a metric via d(x,y) = \\x — y\\. If V is a complete metric space under this norm, then it is 
called a Banach space. The classical Lebesgue spaces are Banach spaces defined by means 
of integral norms. The Lebesgue L p -norm || • || p for 1 < p < oo is given by = ||/||lp : 

= (/ \f(x)\Pdx) 1/p . and the Lebesgue space LP(R n ) is the vector space of all measurable 
functions on R™ with finite L p -norm. (In our case, we consider real functions.) These spaces 
are Banach spaces. The closed unit ball of a Banach space is compact in the norm topology 
if and only if the space is finite dimensional. A Banach space (or even a topological space) 
is separable if it contains a countable dense set. All the classical L p spaces, except for L°°, 
are separable. 

To each Banach space X corresponds a dual space, denoted X*, comprising continuous 
linear functionals. A common notation for the value of the functional A G X* on x G X is 
(A,x). An inequality of the form \(u, f)\ < \\u\\x* \\f\\x expresses the content of continuity 
for a linear functional, and simultaneously serves to define a norm on X*, which turns it 
into a Banach space. The inequality also shows why continuous linear functionals are called 
bounded. 

The duals of the Lebesgue spaces can be identified with Lebesgue spaces themselves. 
For 1 < p < oo, L P (M™)* = L 9 (R 3 ) with l/q + 1/p = 1 and the duality pairing given by 
an integral: (A,u) = j X(x)u(x) dx. For p = 1, the formula gives q = oo. L°° = (L 1 )* 
consists of functions which are essentially bounded, meaning that for u G there is some 
a number ||u||oc M such that \u(x)\ < M off a set of measure zero. H^Hoo, the norm of w, 
is the smallest such M. The spaces LP for I < p < oo are reflexive, meaning they are the 
duals of their duals. L 1 and L°° are not reflexive. 

The Hahn-Banach theorem says that whenever A is a bounded linear functional on some 
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vector subspace V of Banach space X with bound M on V , then there is a (non-unique) 
extension of A to all of X with the same bound. 

For two Banach spaces X and Y , the Banach space XC\Y is the set intersection normed by 
II^II-XTiy — \\ x \\x + II^IIf- The dual space is X* + Y*, normed as ||u|| x*+y* = hrf{|| M i|| x + 
\\u2\\y '■ ui + u 2 = u}. 

Weak topologies. Apart from the norm topologies, there is another class of topologies 
on Banach spaces which appear in this paper. The weak topology on a Banach space X is 
given by neighborhood bases of the origin of the form 

U(Q;{Q,e) = {yeX:\(Z i ,y)\<e: i = l,...,n}, 

for some £1, • ••,£« € X* (a finite set) and e > 0. Neighborhoods of other points in X are 
obtained by translation. Thus, a sequence x±,xi,... converges to x in the weak topology 
if and only if for any £ and e, there is some iV(£, e) such that i > N(£, e) implies that 
Xi G U(x, £, e). Reversing the roles of X and X* gives the weak-* topology on X* . The 
Banach- Alaoglu theorem, which is used a couple of times in this paper, says that a norm- 
bounded, norm-closed subset S of X* is compact in the weak-* topology. If X is separable, 
then the weak-* topology on S is metrizable, so weak-* compactness of S is equivalent to 
sequential weak-* compactness. If X is reflexive, then the weak topology on X is the same 
thing as the weak-* topology when viewing X as X**. 

Appendix B: continuously differentiable convex functions in finite dimensions 

This Appendix contains a brief sketch of the convex analysis results used in the proof of 
Lemma FOl This material can be found in §6.7 of Ref. 44], and also in Ref. 45] with more 
work. It seems advisable to give a proof here because many books on convex analysis do 
not discuss it. 

Here is the general setting in which we work: {/„ : n = 1,2,...} and / are convex 
functions on Br = B(0, R), the open ball of radius R in R™, which are uniformly bounded 
above and below. (5) is the main result. 

1) . The bounds, say a < b, imply that if a; € Br, there is some open neighborhood of x on 
which \f(y) — f(y')\ < c\y — y'\, for c which depends on the neighborhood. This is called 
local Lipschitz continuity. 

For, if x £ Br, some open ball B(x, 25) around x is contained in Br. From the bounds 
and convexity of / it is easy to see that a chord to the graph of f n in B{x, 8) then cannot 
exceed {b — a)/S. 

2) . The pointwise suprememum \J /„, and the pointwise limit lim„ f n (x), if it exists, are 
convex over Br, since lack of convexity can be diagnosed from only three points. From 
these it follows that limsup n f n (x), which always exists, is convex. 

3) . If /„ — > f pointwise, then the convergence is locally uniform: given x and tolerance e, 
there is a neighborhood U of x such that \f n — f\ < e over U , for large enough n. This 
follows from the fact that {f n } and / are equi-Lipschitz. 

4) . Now, if /„ — > f pointwise, then limsup/^(a;; y) < f'(x; y) for any y 6 M n . 

This is essentially a one-dimensional problem, and hinges on the fact that the difference 
quotients in the definition of directional derivative converge monotonically: With 5\f(x;y) : 
= X^ 1 [f(x + Ay) — f(y)] for A > 0, S\f(x;y) I f'{x;y) as A J, 0. To prove the assertion, 
fix e > 0. Clearly, 3Ao > such that A < Ao implies that 8\f{x; y) < f'(x; y) + e/2. Then, 
there is N such that n > N implies that both \f n (x) — f(x)\ and \f n (x + X^y) — f(x + 
X y)\ are less than A e/4, so that 5\ f n (x;y) < 5\ f(x;y) + e/2 < f'(x;y) + e. Due to 
monotonicity of 5\f n (x;y), therefore, n > N e implies that f' n {x;y) < f'(x;y) + e, and 
finally, limsup n fjx\ y) < f'(x;y). 

5) . If f n —> f pointwise, then for x G Br, df' n {x) C df(x) + eB\ for large enough n. If /„ 
and / are Gateaux differentiable, this implies that Df n — > Df pointwise. 
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Since f'(x; y) and f' n {x\ y) are themselves convex functions of y for fixed x, we can apply 
(2) and (4) to see that n > N(y) implies f' n {x\ y) < f'(x; y) + e. But, (3) shows that N{y) 
can be chosen independently of y in the unit sphere. 
6). Applying (5) to f n (x) := f(x — z n ) with z n — > shows that / is C 1 . 



Appendix C: Genericity of X*-nonrepresentable densities 



Proposition C.l. Given (po,vo) £ OF, For any e > and A > 0, there is an X*- 
representable density p with \\p — Po\\x < e and \F[p] — F[po]\ < e, such that \\v[p\p\\\ > A, 
and therefore \\v[p]\\x* > A. 



Proof. Since Thm. 16.21 showed that any A*-representable density is approximable by L c 
representable densities, assume that «o £ L°°. 
Define 



w e := r?(x)sin- = lr)(x) ( ~^" cos — ) > 



(CI) 



where rj(x) is a smooth compactly supported bump function (77 £ C^°(R ) and < r\ < 1). 

The rough idea is that for very large M, if € is small enough, the ground state density 
for potential v + Mwg cannot exploit Mwg because doing so requires too much oscillation 
and therefore costs intrinsic energy with the result that that ground state density is hardly 
different from p$. A slightly indirect approach and possibly a small additional modification 
of the potential is necessary to prove it. The bound (|C2I) below and the Ekeland variational 
principle are key. 

Integrating by parts, 



we(x) p(x)dx 



2£ 



cos ■ 



p71 ^x7 1 dx 



Applying the Cauchy-Schwarz inequality and using that cos 2 x\jl < 1 yields, 



{w t ,p)\ < 2£[ / npdx 



1/2 



riVp^fdx 



1/2 



p\Vr] 1/2 \ 2 dx 



1/2' 



Finally, using inequality ([5]), and overestimating the first and last integrals, 

\(w e ,p)\<£(c 1 F[p]y 2 + c 2 ), 

for some positive constants c\ and ci which depend on total particle number. 

First, assuming I small enough (depending on p ) that ||Mu^/9 ||i > ^ / Povdx =: C3M 



(C2) 



fix M > 1 large enough that (uniformly for small £) 

UMwe + vo'jpoWi > 2A, 



(C3) 



and without loss of generality assume that e < 1/M < 1. A few more upper bounds on £ 
will be imposed. 

From now on, we will be interested only in densities p satisfying the condition 

AS (wi, p) < F[p] + (v + Mw e , p) < F[p Q ] + (v Q + Mw e , p ) = M (w e , p Q ) (C4) 
Define e by 

^ 2 :=e-mm{l,\\v U,\\po\\ x l }<e. (C5) 
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Claim: if I is small enough, then M(wt, po) < e/3, and for any p satisfying (|C4[) . M(wt, p) > 
e/3. The hrst part is simple since (we, po) — > as I — > 0. Given that, the second follows 
by use of the bound (|C2[) . Essentially the problem involves a quadratic in Flp] 1 / 2 . If Ml 
is very small, violation of M(we, p) > e/3 would require such a large F^] 1 / 2 that the upper 
bound in (|C4|) would also be violated. We will not write down explicit criteria, but assume 
i has been so constrained, with the result that F[po] + (vq + Mwe, po) < E[vq + Mwe] + e. 

The ground is thus prepared for an application of the Ekeland variational principle. It 
guarantees existence of a ground-state density p of v : — v + Av + M wg where the pertur- 
bation Av satisfies ||At;||x* < e 1 ^ 2 , and p satisfies \\p — po\\x < e 1 ^ 2 and (|C4[> . Thus, 

F[p] < F[p ] + (v , (po - p)) + (Mwi, {po - p)) + (Av, ( Po ~ p)) 
< F[po\ + e + I + I 1 ' 2 < F[po\ + 3e. 

This shows that we get p close to po in L 1 norm and intrinsic energy as claimed in the 
statement of the Theorem. The upper bound on F[p] suffices due to lower semicontinuity. 

All that remains to be checked are the sizes of \\vp\\i and Calculate using (|C3[) 

and (|C5|) : 

||(u + Av + Mwg)p\\i > \\(vo + Mw e )po\\i - ||(^o + Mw e )(p - p )||i - ||A«p||i 

> 2A - (M + |K||oo)e 1/2 - ||A«||x' ||po||x - \\Av\\ X '\\p - Po\\x 

> 2A - (1 + e) - e - e > 2A - 4. 

Finally, notice that |M|x* > \\ v p\\ i/\\p\\x and that \\p\\x > 1- □ 

This proposition will be used here as a springboard toward the conclusion that failure 
of A*-representability is topologically generic in a nontrivial sense. Topological genericity 
is meant here in the sense of the Baire Category Theorem. Recall that in a complete 
metric space S, a subset is generic if it is a countable intersection of open dense sets. A 
subset is meager if it is a countable union of nowhere dense (empty interior) closed sets. 
The complement of a generic set is meager. Obviously, a countable intersection of generic 
subsets is itself generic. Though this notion of 'generic' is not without flaws, it is commonly 
accepted as a criterion of just what the name implies. At any rate, a generic set is dense. 
To apply the concept, we need a reasonable metric which makes Jn a generic set in its 
completion. Thus, equip Jx with the metric 

d(p,p') = \\p-p'\\x + \F[p]-F[p'}\. (C6) 

This metric is actually just the graph norm of F as graphic clxt an observation which 
is key to identifying the completion of Jn with respect to d with epi F C X x R. Pictorially, 

(Jx,d) — (Jn,<I) 



graph F > epiF 

The identification is established by showing that epi F is a closed subset of X x M (hence 
a d-complete space), and that JTjv is dense in epiF. The first is a simple consequence of 
lower semicontinuity of F. The second of the fact that F is unbounded above on every 
A-neighborhood of any point in Jx'- take (p, z) £ epiF, z > F[p], and e > 0. Then, there is 
p' satisfying \\p' — p\\x < e and F[p'] > z. Over the line segment [p, p'], F goes from F[p] < z 
to F[p'] > z, hence must pass through z. This shows that there is a point (p" , z) £ graph F 
within d-distance e of (p,z). Hence, graph F is d-dense in epiF. In fact, Jx is actually 
generic in Jx ■ For, since epiF is closed, so is epi(F + a) for a > 0, since this is just a 
translate of epiF inlxl. The argument in the previous paragraph shows that epi (F + a), 
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as a subset of (epiF,d), has empty interior. Therefore U n >iepi(-F + * s meager in J^. 
But this last is all of epiF except for graph F. This result is crucial to the endeavor. A 
property P defined on Jn is generic if and only if P-and-in-j7}v is generic, so that the status 
of elements in Jn \ Jn is irrelevant. 

Paraphrased in terms of d, the Proposition says that neither ||u[p]p||i nor (with- 
out loss set equal to +00 outside the effective domain) is upper semicontinuous with respect 
to d. Denote by A a the set A a : = {p E Jm : < ct], of densities with representing 

potentials having A* norms not exceeding a. According to Section A a is closed with re- 
spect to A norm, a fortiori with respect to d, due to lower semicontinuity of ||i>[']||-X"* ■ On the 
other hand, the Proposition shows that A a is nowhere dense. Consequently, 1Z = W^ =1 A a 
is meager, being a countable union of closed nowhere dense sets. Since 1Z contains every 
A*-representable density, not being AT*-representable is generic. 
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